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摘 要 : 硬 峭 值 追踪 算法 是 压缩 感知 中 一 种 重要 的 重 构 算 法 ， 该 算法 的 本 质 是 一 个 最 小 二 乘 问题 ， 存 在 复杂 度 高 ， 
收敛 性 差 、 运 行 时 间 长 等 缺点 。 引 入 Nesterov 方法 用 以 优化 稀世 解 的 凸 松弛 现象 ， 引 入 逐次 松弛 迭代 法 优化 传统 硬 
阅 值 线性 方程 组 ， 理 论证 明 优 化 结果 具有 良好 的 收敛 性 ,仿真 实验 说 明 优 化 后 的 算法 有 效 降低 算法 复杂 度 ， 降 低 了 
运行 时 间 。 
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Research on optimization hard threshold tracking algorithm 


Chen Xuan!, Pan Chun Pin!, Long Dan? 
(1. Zhejiang Industry Polytechnic College, Shaoxing Zhejiang 312000, China; 2. Zhejiang University, Hangzhou 310058, 
China) 


Abstract: The hard threshold tracking algorithm is an important reconstruction algorithm in compressed sensing, which is 
essentially a least squares problem and has the disadvantages of high complexity, poor convergence and long running time. 
This paper introduced the Nesterov method to optimize the convex relaxation phenomenon of sparse solutions, and 
introduced the successive relaxation iterative method to optimize the traditional hard threshold linear equations. The theory 
proves that the optimization results have good convergence. Simulation experiments show that the optimized algorithm 
effectively reduces the complexity of the algorithm and the running time. 
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0 引言 司 的 角度 其 进行 了 研究 ,取得 了 一 些 成 果 。 文 献 [14] 提 出 S- 闭 
一 包 的 分 块 相干 性 和 约束 等 距 特 性 (RIP) 概 念 ;根据 元 余 字典 分 
压缩 感知 理论 中 突破 传统 的 Nyquist 定理 的 严格 限制 ， 岂 相 干 性 的 特点 ,对 更 新 支撑 集 部 分 对 算法 进行 改进 ,仿真 


实现 了 对 信号 的 采样 能 够 在 较 低 的 采样 率 下 进行 ,并 同时 完 验 说 明 有 效 的 降低 了 算法 的 收敛 性 ;文献 [15] 提 出 一 种 基于 
成 对 信和 号 的 压缩 ,在 重 构 一 端 通过 非 线性 化 的 方式 进行 信号 SHTP 算法 的 联合 重 构 算法 来 求解 分 布 式 压缩 感知 问题 , 仿 
的 重 构 。 压 缩 感知 理论 的 提出 使 得 信号 在 采样 过 程 中 不 再 主 结果 表明 该 算法 具有 较 大 的 信号 重 构 噪 声 比 和 较 小 的 平均 支 
要 依赖 信号 的 宽度 ,这 样 能 够 有 效 的 降低 了 信和 号 传输 和 存储 撑 势 误差 ,可 实现 信号 值 的 精确 重 构 ; 文 献 [16] 提 出 一 种 压缩 
过 程 中 所 需要 的 代价 。 因 此 研究 如 何 能 够 高 效 的 恢复 原始 信 采样 硬 闵 值 的 重 构 算 法 ,该 算法 在 去 除 原子 时 结合 人 硬 阔 值 思 
号 成 为 了 重 构 算法 的 主要 方向 ,目前 在 压缩 感知 中 的 重 构 算 想 ,保证 了 每 次 迭代 更 加 的 精确 ,仿真 实验 说 明 算 法 具有 更 好 
法 有 则 松弛 算法 , 贪 禁 追踪 算法 和 组 合 优 化 算法 。 凸 松弛 算法 9 重 构 性 和 抗 躁 能 力 ;文献 [17] 提 出 了 一 种 新 的 部 分 硬 阔 值 
观测 次 数 少 ,复杂 度 偏 高 的 特点 ,其 代表 有 基 追 踪 算 法 [9 (PHT) 算 子 , 它 类 似 于 硬 闵 值 算 子 ,限制 了 支持 集 在 一 次 迭代 
(basic pursuit, BP), BB BE Ht EG Fa Wu E Œ $E A Ulgradient 中 可 以 改变 的 量 ,仿真 实验 说 明 该 算法 提高 算法 效率 ;文献 [18] 
projection for sparse，GPSR), 最 小 角度 回归 法 [Sl(least angle 提出 了 一 种 扩展 的 稀疏 恢复 理论 并 通过 在 硬 闵 值 算 法 中 使 用 
regreesion，LARS); 贪 焚 追 踪 算 法 具有 算法 具有 节点 ,速度 快 ， 取得 了 较 好 的 效果 ,能 够 有 效 的 恢复 信号 ;文献 [19] 提 出 了 一 
复杂 度 低 的 特点 ,其 代表 有 正 交 匹配 追踪 算法 外 (orthogonal 种 复杂 的 洛 伦 效 迭代 硬 阔 值 算 法 ,仿真 实验 说 明 该 算法 能 够 
matching pursuit, OMP ), 正 则 化 正 交 匹配 算法 ll(regularized ” 有效 的 恢复 信号 ,具有 较 好 的 效果 ;文献 [20] 将 硬 闪 值 算法 用 
orthogonal matching pursuit,ROMP), 分 段 正 交 匹配 追踪 算法 于 从 少量 线性 测量 中 最 优 地 恢复 稀疏 信号 ,然后 对 CS 进行 去 
[ll(stagewise orthogonal matching pursuit, STOMP) 等 ;组 合 优 品 和 传统 去 品 处 理 的 比较 研究 ,该 方法 能 够 有 效 的 降低 噪声 
化 算法 主要 是 稀疏 傅 里 叶 变 换算 法 02(sparse Fourier 对 信号 的 影响 ;文献 [21] 提 出 在 硬 阔 值 算 法 中 将 黎 曼 优化 算法 
transform, SFFT). RME p] ftj iterative hard thresholding H F HK EAE RE PK S RI IE AE , Deu FE gu p EE SE £3 PE CSS P AE 
IHT) 是 一 种 凸 松 弛 算法 ,优点 是 适合 大 规模 高 维 数据 人 处理 和 和 础 秩 和 矩阵 。 经 验 评估 表明 ,算法 能 够 从 几乎 所 需 的 最 小 测量 
实时 数据 处 理 ,缺点 是 测量 和 矩阵 在 有 限 等 距 特性 条 件 下 需要 数量 中 恢复 低 秩 矩阵 ;文献 [22] 提 出 了 一 种 基于 广义 期 望 最 大 
满足 最 大 奇异 值 小 于 1 的 限制 条 件 ,因此 导致 了 算法 对 对 测 化 (GEM) 噪 声 稀 玻 重 构 的 降 噪 硬 阔 值 方法 ,取得 了 较 好 的 效 
量 和 矩阵 的 尺度 要 求 高 ,运算 运行 时 间 长 ,复杂 度 高 。 学 者 从 不 果 。 
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Ak x TE OC BR [23] 2& F 88 I] [t C25 9T A HS A uli E, 5| 
Nesterov 和 逐次 超 松弛 替代 法 来 进行 改进 ,并 证 明了 改进 后 的 
和 法 具有 良好 的 收敛 性 ,在 仿真 实验 的 指标 测评 中 取得 了 较 
好 的 效果 。 


1 AREER? 


所 有 的 压缩 重 构 算法 都 必须 遵守 限制 自 同 构 性 质 
(RIP), ik Aem 的 s 阶 的 自 同 构 系数 5.=5.(4) 定义 为 满足 
式 (1) 的 最 小 的 5。 

4-8) xls] Ax Is a+ I x lk (1) 

对 于 任意 的 5 稀疏 的 向 量 xeR"»。 文献 [20,21] 证 明了 当 A 
满足 了 RIP 条件 后 ,这 些 算法 都 是 收敛 的 。 

硬 阔 值 追踪 算法 结合 了 IHT 的 硬 冰 值 策略 和 CoSaMP 的 
投影 算法 ,迭代 公式 如 下 : 

Sn 2 L Go, A (y - Ax,)) Q) 


Xan -argmin || y— Az |. supp) c Sn 3) 


式 (3) 中 的 supp(x") RIR o^ BSc PESE. LOO 表示 x 绝对 
值 中 最 大 的 * 个 元 素 的 支撑 。 硬 阔 值 的 欠 代 步 又 ( 即 式 (3)) 的 
本 质 形式 是 一 个 最 小 二 乘 问题 。 随 着 压缩 感知 中 的 信号 重 构 
的 问题 规模 逐渐 增 大 ,采用 直接 解决 最 小 二 乘 的 问题 就 会 非 
常 的 耗 时 ,因此 ,文献 20] 中 提出 的 一 种 快速 硬 阔 值 的 算法 能 
够 有 效 的 解决 大 规模 下 的 稀疏 恢复 的 问题 , 它 采 用 的 梯度 法 
来 解决 先 代 中 的 最 小 二 乘 问题 , 当 使 用 梯度 法 的 时 候 , 取 步 长 
为 1, 得 到 : 


IERI x, -x]h (4) 
其 中 :x JETER s PARIE 4x=y 的 向 量 , Udo. 是 快速 硬 
闵 值 的 算法 产生 的 迭代 序列 。 从 上 面 的 研究 中 发 现 , 式 (3) 一 方 
而 是 一 个 待 优 化 问题 ,采用 优化 效果 更 好 的 Nesterov 可 以 提 
高 优化 的 效果 ,从 另 一 个 方面 来 看 式 (3) 也 是 一 个 正定 的 线性 
方程 组 ,需要 使 用 线性 解决 办 法 更 好 的 逐次 松弛 替代 法 解决 
问题 ,本 文 阐述 并 证 明了 式 (3) 通 过 以 上 两 种 算法 的 改进 得 到 
的 算法 仍然 是 收敛 的 。 
2 BOB E) EDO ER RE 
2.1 Nesterov 方法 

该 方法 是 一 种 用 以 求解 加 速 求解 光滑 优化 问题 的 算法 ， 


IA 


Ixa xl 


min fœ (5) 
其 中 : f 是 Lipschitz f EE L7 LC , 即 满足 
l| VF Ga) - Vf Go) |b S L| y 7x3 lbs xx» € domCf) (6) 


算法 步 又 如 下 所 示 。 
算法 1: Nesterov 算法 
1: 初始 化 : y =x, =l 
2: while "not converage" ,do 


xe e yt -AVO 

Aa - 50-6048 84-8) 
yet oae f Gt cat) 
8,756 8 74-8) 


k=k+1 
endwhile 
最 终 的 迭代 结果 x* 
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算法 中 的 及 是 步 长 参数 , 设 定 =, 文献 [12] 证 明 存在 如 
下 的 关系 


2L|| y? — y |È 
(k 1? C) 


F -fO 


其 中 的 六 就 是 问题 (3) 的 解 且 满足 ESO 。 当 式 (3) 使 用 算法 
FER CA f CO | y -Ax|B , A, AG. S) e R^ , 根 
dis A 具有 限制 自 同 构 性 质 得 到 : 


lA xls a-5)| x|B«2. Yx e Rs (8) 
KEAR PRI M Ax B2 LB EA db 2 ,因此 有 
I| VF G3) - V£ G3) lb SE AA, C08 722 Il 2 IE — x Dl (9) 


2.2. 逐次 超 松弛 和 返 代 法 
式 (3) 其 本 质 也 是 求解 下 列 线性 方程 
AAXTAD (10) 
A, 2 AG,S,) ER" ,因为 4 具有 限制 自 同 构 性 质 , 所 以 , AA, 
是 正定 矩阵 ,把 矩阵 AA, 分 裂 成 人 4 -D-(D- AA) -D-L-E , 
其 中 , D=diag(4d1,d,,…d,) ,4>0 是 对 称 正 定 矩 阵 44 的 对 角 部 
分 ,-L 是 和 4 的 严格 下 三 角 部 分 。 令 AASM-N 


M =L(D-_wL) 
Ww 


N=L(0-wWD+wE),w=0 
w 


则 得 到 逐次 超 松 弛 迭代 法 
a9 2 wD^ (Lx? + E79 + A*y) + wx (11) 
即 
x 2(D-WL)*(02w)DWE)x? -w(D-WL)'Ay (12) 
其 中 : w AERA. E iE ICE X 02) 8935 (OB FE 


G, - (D- wL) (0-w)D - wE) ,上 述 迭 代 需 要 显示 计算 44 , 若 在 
58 k 2p, ZO eaten, rO -b— Axfe» , 则 可 得 不 需要 显示 计算 44 。 
算法 步骤 如 下 所 示 。 

算法 2 逐次 超 松 弛 迭代 法 算法 

1: ZO =X ，70 = Ay — A. x07» 

2: Zi» zZ0 +ő;e; 


rid 2r 一 51aj 
ô; 2 wa;ri?[d;, j sl.-n 
x0 = Zo 
3: 输 出 迭代 结果 xt 
2.3 ”本文 算法 的 收敛 性 分 析 
本 文 算 法 步骤 如 下 所 示 。 
算法 3 本 文 算法 
1: 初 始 化 : x 
2: while "not converage" ,do 
Spa =L G, c A'(y- Ax.) 
Xna = QS t A (Y = Ax), 


xà -(0D-WL)(1-W)D-c-WD)xGi a 

*w(0-WL)? A; y " 

k-k4l 

iflxa-xl«e 转 步 又 3, 和 否则 转 步 又 2 
endwhile 

3: 输 出 结果 x 
下 面 的 过 程 就 来 证 明 算 法 3 是 具有 收敛 性 的 
定理 1 i AER“, xeR,yeR" , Aze A BACTUTE EIERE, 
则 当 0«w «2 的 时 候 ,迭代 方法 3 是 收敛 的 。 
证 明 设 14 是 逐次 超 松弛 迭代 法 方法 的 友 代 入 


ru 

Tr 
A 
到 


i 
4g 5 
B 


任意 一 个 特征 值 ,XsRs 为 其 相对 应 的 特征 向 量 , 则 有 等 


(0—w)D- wE)x 2 ACD -WD)x ,因此 
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-Wx Lx») ,因为 AA, 对 称 正定 


(1-W)x' Dx -Wx'Ex = A(x" Dx 
所 以 xDx2420 


x'Lx- a ^ fli 
x'Dx-(x'Lx) =a- fli 
XAAx-x(D-L-D)x-q-2a 


EHH, æ, BeR 


HB adc w+ we wpi 
TE q—wa —wfi 


wl aP +w? p? 
可 得 | 外 [g Eu EN P 
当 0<w<2 的 时 候 , 有 
[q-wq-a)P -(q- wa} » qw2.— w)2a —-q) «0 [A] JG, | A]«1 
即 有 Pp(Gw)<1 
因此 ,和 欠 代 方法 是 收敛 的 
因此 接 下 来 证 明 算 法 3 具有 的 收敛 性 
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引 理 1 根据 文献 [20] 中 的 假设 4 
U 是 指标 集 。 因 此 对 于 任意 的 上 铅 量 v 来 说 ,存在 
I| (4 — A4), so ||v |l; 


mt 
o 


HH, JUUsuppQ) et 


FRW AE supp) c S, In] 5c z LEE s PSI EAE Y 
量 *, 有 


(x= x2) < z || xx; lell z ll 


F Ax 是 {hz,supp(z)c5,} 上 的 yi 


证 明 
有 


近似 最 佳 ， 


(y- Ax;, Az) Oif. supp(z) c S, 
因为 >=4x, 式 (15) 重 写 为 

(Ax-x,,Az) - 0,if. supp(z) c S, 
将 式 (16) 和 引 理 1 RARAS), EX U, 9 supp(x x, ) 
(x-x,.Q 一刀 A)z) 


x; lell z ll 


(x-x,2)- 


(x-x.[Q -AA z], ) S s | x 


引 理 3 在 算法 3 Pos FAE Fa a s BW AE y= Ax 的 向 量 


x Uses 满足 
lO =s, lb 5, | ox — x lb 
K X, + A'(y — Ax, as 
证 明 a E ， 
在 两 边 同时 减 去 上 % +A O- Ax Ds, s.s 有 ,得 到 
l| Gc, +A" Qr Ax, Js, s... lb 
«| — A AX, — 2), s... ll 
2| G- x5, llb Il — A AG, -ass 因此 ,得 到 
lx), bsll (X — A DG, 735, s... l 


WE S = supp) ,得 到 


5382 在 算法 3 中 ,定义 lim 六 = 总 。 因 此 对 于 任意 的 * 


(15) 


(16) 


hn 


(17) 


引 理 4 设 雹 是 由 逐次 交代 方法 产生 的 序列 , 则 
iae: i o 

le-a es ge] 
定理 2 设 定 算法 3 中 产生 的 序列 a 满足 


l| xi — EE 7 E x0 — l, AE A TF EE E s PRU AL AE. y Ax 
量 x, 因 此 算法 3 PRIE x V 


lx, -xlbs e" Ilx -xl 


im] 
pii 


1 
p=(1+6,)(6,, +70) UE 


+y (D) * Ó;, 
证 明 本 文 可 以 得 到 
Il [tom 一 20s ， l= (Xa =X, Qa —X)s,, j= 


(Xa E T E X9) t (x X, Qa X55) 


的 向 


(18) 
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注意 到 w=, ,因此 ， 


X Xn = Xs, < YO a 7 Xa db xl Ga 7 x05, lb 
通过 引 理 4 
Xa 7X Q8 7X)s, < Cs xa =x ll xl Goa = xs... ll 
化 简 得 至 


Il Qua —-X)s,. 1 ll; 
€à, | xa —xclb YO xta xa lb 
X (&, X || xia =le 0l xa x lle 


根据 4 的 限制 自 同 构 i ibd 
5, l^ "L1 


(19) 


| xs xls 


, l|(D-WL)^ (A -W)D - WE )xki 


n+l 
+w(0-WL)' Ay- y | a sl X)s,, 一 站 


1+0， 
1-6, 


|| (D- WL) (Q -W)D+WE )xzi 


Ità, 
+w(0-WL)' A y-x)s, lB + 


Osa lE 


g 1+, 
l- ô, 


"E 


I 


| Gs XIE + lie Xu —X)s,, IE 


到 
I| xta =x BS G5, + A Qr 7 Ax s. —xlB 
S| Œn x A O= Ax, DB s.s IB SIE 7 A677 ns, ME Ds, BB 
< ó$ [lx =x Mb +l Gua =s, lÈ 
两 边 开 方 得 到 
| — xl =s, ll 
Ac EA RRCA 9) 65:581 


Il eom —X)s,, ll; 


<[(6, - y (D)) p 5 -Ó,Q0]l x, —xlb + 
(8, ero s MAD Ds lk 


由 引 理 3 得 到 将 式 (18) 代 入 上 式 中 得 到 


| (Xn —X)s,, ll 


€à, | x- x, lb +I Ga 


1-ó;, 
1-ó, 


2s 


€ [0 eXe +70) 


*óy()] | XQ —X Il: 


最 后 可 以 得 到 
=x G6 =s, lb + 
| Ga =s, lbs ell x, 7x 
忆 此 得 到 : Mels el xl, 因此 本 文 算法 证 明成 立 。 


3 ”实验 仿真 


3.1 算法 重 构 性 能 对 比 
为 了 进 一 ga qam eu da 
256, RREA K 的 随机 信号 *, 设 定 @ 73 MxN 为 随机 高 斯 
阵 , 测 量 值 为 yY=@x ,实验 信号 的 稀疏 度 比 值 设 定 为 [0.05,0. 
之 间 , 重 构 精 度 以 误差 率 fli HII 作为 判定 比较 标准 
设 定 当 误差 率 小 于 1e-004 的 时 候 , 说 明 信 和 号 重 构成 功 ， 当 误差 
率 大 于 等 于 1, 则 说 明 信 号 重 构 失 败 ,将 OMP,SP,CoSamp,IHT 
与 本 文 算法 在 重 构 性 能 和 运算 复杂 度 两 个 指标 上 进行 对 比 。 
设 定 实验 运行 300 次 ,对 这 5 种 算法 的 实验 数据 结果 按照 取 平 
均值 进行 计算 。 图 1 和 2 显示 了 这 5 种 算法 在 不 同 的 稀疏 度 
下 的 重 构 精 度 以 及 算法 所 需要 的 时 间 。 
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10 poer Y ——3À r : [100,200], £47 500 次 实验 后 的 平均 数据 。 表 2 显示 了 5 种 算 
os] M | 法 的 重 构图 像 的 PSNR 值 ,从 表 2 中 可 以 发 现 ,本 文 算法 的 重 
o.8 NS -了 ] 构图 像 的 PSNR 值 相 比 CoSamp 和 IHT 提高 了 0.2-0.6dB, 比 
MN x SS | OMP 和 SP 提高 了 0.8-0.9 dB. 
= 95] Ve NR. ] 表 2 不 同 算法 下 的 重 构 性 能 和 时 间 对 比 
ka 0.54 v A. ^ ps ] Table 2 Comparison of Reconstructing Performance and Time under 
ii +, 
E o P IW X ] Different Algorithms 
dmt JN [n] E 
& 034 NSW. M-100 M-200 
0.24 Te \E -< 算法 
. A V V PSNR 时 间 PSNR WA 
了 ES w ] OMP 28.94 6.99 31.07 9.76 
0.0 - T TA A ecc 
0.05 0.10 0.15 0.20 0.25 0.30 SP 29.13 7.13 31.23 21.58 
Mii SE CoSaMP 29.71 16.52 31.29 42.95 
图 1 不 同 算法 下 的 重 构 精 确 率 对 比 IHT 29.82 4.48 31.24 4.57 
Fig. 1 Comparison of the Reconstructing Accuracy under Different 本 文 算法 29.98 16.75 31.86 44.07 
Algorithms 为 了 测试 该 算法 在 确定 型 观测 和 矩阵 下 精 够 重 构图 像 , 观 
HL | ; 测 矩阵 部 分 选取 部 分 哈达 玛 矩 阵 ,实验 结果 如 图 3 所 示 。 
0.94 TZ omP A 36 T T T T 
—€-— CoSamp m —m— OMP 
985 ODE 84] |-*— SP 一 本 
0.74 SUC | J] | —— CoSamp T d 
AX ee ad mq, —A4 
E 064 yU .ct ] —— 本 文 算法 AN 
t£ EA e-e 一 
ES -—v e-e 
T d QUI i ] S =] e A 
EN 0.44 atvett J z "i E | 
$ o3 NU ] 9 Jj 
aa] . n7] 26 jet 
A ve m-i E 
M d oro-o-o-or-oroooo nme nt = 24 
0.05 0.10 0.15 0.20 0.25 0.30 as w xm i x 
sU 预测 数目 
图 2 不 同 算法 下 的 平均 运行 时 间 图 3 不 同 算法 下 的 哈达 玛 矩 阵 下 的 PSNR 值 对 比 
Fig.2 Average Operation Time under Different Algorithms Fig.3 Comparison of PSNR Values under Hadamard Matrix under 
图 1 显示 了 在 不 同 的 稀疏 度数 值 下 的 5 种 算法 重 构 精 度 Different Algorithms 
的 对 比 ,从 图 中 发 现 随 着 稀 琉 度数 值 不 断 增 大 ,5 种 算法 的 重 图 3 中 显示 了 随 着 观测 数目 不 断 增加 ,五 种 算法 的 重 构 


构 精度 都 在 不 断 的 减少 ,从 图 中 的 曲线 中 发 现 本 文 算法 是 在 。 图像 的 PSNR 数值 也 不 断 提高 ,从 整体 上 看 本 文 算法 均 好 于 其 
WEB 0.2 位 置 才 开 始 下 降 ,是 所 有 曲线 下 降 最 晚 的 ， 也 几 种 算法 ,采用 部 分 哈达 玛 矩 阵 也 能 精确 地 重 构图 像 。 因 此 
这 说 明 本 文 算法 在 重 构 精 度 方面 具有 一 定 优势 。 图 2 显示 了 ” 本 文 算法 能 够 有 效 地 重 构 信和 号。 

各 种 算法 平均 运行 时 间 , 伴 随 着 K 的 不 断 增 大 ,本 文 算法 的 运 
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行 时 间 相 比 于 其 他 算法 是 最 短 的 。 zl 
3.2 算法 抗 躁 性 能 对 比 重 构 算法 一 直 都 是 压缩 感知 中 的 重要 研究 方向 ,本 文 针 


为 了 进一步 说 明 该 算法 具有 的 抗 躁 效 果 , 本 文选 取 了 选 ”对 硬 阔 值 追踪 算法 自身 存在 的 不 足 ,提出 了 引入 Nesterov 和 
Bp E k 值 为 30 的 上 述 信号 , 设 定 。 为 干扰 噪声 ,观测 向 量 — 逐次 超 松 弛 进行 改进 ， 仿 真实 验 说 明 本 文 算法 能 够 有 效 的 提 
表达 式 为 Y=qx+e ,干扰 噪声 设 定 分 别 为 0.01,0.05,0.1,0.15, 分 。 高 收敛 性 ,降低 运行 时 间 , 提 高 了 重 构 效果 。 
别 对 OMP,SP,CoSamp,IHT 以 及 本 文 算法 进行 抗 噪 性 对 比 , 实 
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